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Abstract 

We will review some of the theoretical progresses that have been recently done 
in the study of slow dynamics of glassy systems: the general techniques used for 
Tij- ■ studying the dynamics in the mean field approximation and the emergence of a pure 

. dynamical transition in some of these systems. We show how the results obtained for 

I a random Hamiltonian may be also applied to a given Hamiltonian. These two results 

open the way to a better understanding of the glassy transition in real systems. 

■ 1 Introduction 

a 

^ I Many systems (among them glasses, rubber, spin glasses...) show at low temperature a 
5^ ■ very slow approach to equilibrium. Our aim is to understand this kind of behaviour, in 
O ■ particular the peculiar properties of the glass transition. 

At the present moment in the framework of the mean field approximation the static of 
spin glasses is well understood. We start to have a reasonable understanding of the dynam- 
^ I ics. We hope that these progresses will lead to a better understanding of the behaviour of 
c3 ; glasses. However this goal seemed to be unreachable for two reasons. 

• (a) The behaviour of spin glasses and real glasses is rather different near the transi- 
tion. 

In real glasses a very interesting phenomenon happens: if the temperature is de- 
creased fast enough, the system goes into an amorphous state that has an extremely 
large mean life; moreover under very slow cooling some of these systems go into an 
ordered crystal phase. The crystal phase does not exist in general, it exists only 
if the parameters of the interactions among the atoms are intentionally chosen in 
such a way that this phase is energetically favoured. Independently from the possible 
existence of the crystal phase, the dynamics in the glassy phase becomes extremely 
slow at low temperature. 

Moreover the internally energy in the amorphous phase depends on the cooling speed 
in a strong way, also for very slow cooling. These effects have never observed in spin 
glasses. The spin glass transition can also be characterised by the divergence of a 
static quantity (the non linear susceptibility) while no anomalous behaviour has been 
observed in the static quantities in real glasses. 
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• (b) From the theoretical point of view the two systems seem to be extremely different. 

Spin glasses are systems with a random Hamiltonian. They are studied by aver- 
aging over the disorder in the Hamiltonian using ad hoc methods, e.g. the replica 
technique. Real glasses have a given Hamiltonian and the disorder is not present in 
the Hamiltonian, but comes out from the freezing of some of the relevant degrees of 
freedom. 

Recently there have been some progresses which show that the previous difficulties are 
less serious than it was believed. 

• (a) It has been found that in some disordered systems (e.g. modified spin glasses) 
there is a dynamical transition which is characterised by the divergence of the time 
needed for equilibrate the system. Near this transition temperature the static quanti- 
ties do show no anomaly. The temperature at which the equilibration time diverges is 
higher than the temperature at which a transition is present for the static quantities. 
This static transition cannot be observed because the system takes to much time to 
equilibrate. 

• (b) In some cases the replica method can be applied also to system without intrinsic 
disorder in the Hamiltonian. The behaviour of those system is very similar to that of 
really disordered systems, apart from the possible existence of a low energy ordered 
crystal phase 

Some of these effects cannot be understood using only the tools of standard equilibrium 
statistical mechanics, because they are non- equilibrium phenomena. In this note I will 
present a coherent picture of the dynamics of these systems. The results have been obtained 
in these recent years mainly in the framework of spin glasses theory [jl], |2[ and I will show 
how they can be extended to other systems like glasses. Some of the statements are well 
proven, while other are still conjectural. 

In the second section of this note I will give a general qualitative description of the 
dynamics based on the hypothesis that systems evolve in time jumping from one to an 
other quasi-equilibrium state. In the third section I will review some very interesting 
results obtained by a direct computation of the dynamical evolution of the system using 
more powerful methods; I will also compare these results with those obtained in the previous 
section. Finally in the fourth section I will address to the old problem of comparing the 
behaviour of a system with random Hamiltonian (like spin glasses), with that of system 
with a fixed Hamiltonian (like glasses). This comparison will be done in a model system: 
one finds that the properties of the system with fixed Hamiltonian are very similar to that 
of the systems with random Hamiltonian. The only difference is the possible existence of 
a crystal phase for specific choices of the parameters of the fixed Hamiltonian; the crystal 
phase does not exists for the random Hamiltonian. Some brief conclusions are presented 
at the end. 



2 Local Equilibrium States 

We consider Ising spin models in which there are variables cxj, which take the values ±1 
Q There are many possible kinds of Hamiltonian that we can write down. In the simplest 
case the interaction involves only two spins, 

^Similar consideration can also be done for interfaces or manifolds in a random medium Istl. 
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Hj{a) = - XI Ji,k(^i(^k, (1) 

i,k 

and all the pairs i, k are equivalent. 

If the J are randomly distributed Gaussian variables, with variance we obtain the 
SK model, otherwise one has a different model. 

At low temperature these models have a corrugated free energy landscape, with many 
local minima, separated by high barriers. In this situation general arguments imply a very 
slow dynamics because the system may be trapped in a valley and it takes quite a long 
time to escape from it. 

The jumping from a valley to an another valley is a controlled by the height of the 
barriers and the time needed (neglecting prefactors) is 

r = exp(/3AF), (2) 

where AF is the minimum barrier in free energy that the system has to cross in going from 
one valley to an other valley 0. 

The best characterisation of a valley is a region of the phase space in which the 
system spends a long time. In this case it reasonable to define the local magnetizations in 
the valley a as 

= {oi)^, (3) 

where the average is taken inside the valley a. 

In long range models these magnetizations satisfies the mean field equations. Neglecting 
the Onsager reaction field 0, they are are 



rrii 



th{(3Y.Ji,kmk) (4) 



Alternatively we can define a free energy as function of the mj: 

F[m] = -J2AkiT^i'^k-Tj2simi), (5) 

i,k i 

where the local entropy is simply given by 

1+m 1 + m 1-m ,1 - m 
s{m) = — ln(^— ) — ln(^— ). (6) 

The local minima of the free energy F[m] are solutions of the mean field equations . 

The free energy landscape can be characterised by the structure of the set of the so- 
lutions (which we will label by Greek indices). Generally speaking the most important 
parameters are the free energy of the solution (/", the local magnetizations mf in a given 
solution, the overlap among two solutions a and 7 (g^^^ = Y.i=i,N ^^i'^nj) and the self 
overlap {qEA = (la,a), which in most of the models is independent from the solution. 

At equilibrium it is reasonable to assume that different valleys may be populated, and 
the probability that the system is in one of this valley is given by 

Wa oc w{fa) = exp(-/?/„). (7) 



^This is not the only method for having a slow dynamics. For alternative possibilities see [s). 
•^Sometimes one uses the terminology local equilibrium state or quasi-equilibrium state. 



If we add the Onsager reaction field we obtain the TAP equations. 
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It is evident that 

E^" = l- (8) 

a 

In many disordered systems the number of valleys as function of the free energy (A/'(/)) 
increases as 

Afif) ~ exp(y(/ - /o) + 0((^^), (9) 

in the region where 

l«/-/o«V^ (10) 

where V is the volume of the system. 

At a given temperature there are two possibilities 

• (a) If y < /3 the integral 

I dfU{f)w{f) (11) 

is dominated by / near /q. Only few valleys dominate the sum in equation (H), 
although an infinite number of them give a non zero contribution. In this situation 
we say that the replica symmetry is broken. 

• (b) If y > /9 the integral 

j dfU{f)w{f) (12) 

would be divergent if we neglect terms of 0{{ ^~y^ ). In reality the integral is dom- 
inated by the region where / — /o of order of the size V of the system. In this case 
the number of valleys which dominates the sum in equation (|[) is exponentially large 
and each of the valley has a weight which is exponentially small. In this case the 
replica symmetry is not broken. 

Let us consider the case where the valleys are separated by very high barriers (e.g. 
diverging with A^) in the region where the temperature T is smaller than T^. Depending 
on the nature of the problem we may enter or in the region (a) or (b) when we decrease 
the temperature from above to below Tq. 

If we enter in the region (a), as it happens in the usual SK model for spin glasses, a 
phase transition is present from the equilibrium point of view at T^. 

On the contrary, as happens in other spin glass models if we enter in the region 

(b), no phase transition is present from the equilibrium point of view at Tp, and a static 
transition is present only at smaller temperature, Tr, where y becomes smaller that (3 
and we pass from region (a) to region (b). This transition can be easily understood. 
At temperatures greater that Tr the valleys populated at equilibrium have a free energy 
density greater than minimal possible because of entropic effects. The relevance of these 
entropic effects disappear by decreasing the temperature and for temperatures less than 
Tr the valleys have the smallest possible free energy. 

In this last case we have two transition one for the statics and the other for the dynamics. 
If we quench the system at temperature smaller than Td coming from an high temperature 
region, the internal energy (for an infinite system) does not go to the equilibrium value 
and the system remains trapped in a metastable state 

In this case for the infinite system independently from the speed of cooling we always 
find the energy of the metastable states and therefore we do not observe any strong depen- 
dence of the energy on the cooling. This is an artefact of the mean field approximation, 
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which is correct only for infinite range forces. We shall see later how this behaviour may 
be changed for more realistic model. 

It is also possible that there is an other isolated solution to the mean field equation, 
with a free energy density Fc smaller than Fq and there are no solutions in the region 

Fc<F < Fo. (13) 

In many cases this isolated solution describes an highly ordered state, which we call the 
crystalline state. The existence of this states does not change the properties of the system 
in the region where f > fo- 

From this point of view in order to compute the approach to equilibrium one should 
evaluate the free energy barriers which separate one valley from an other valley. This 
computation is rather difficult, especially if we take care that the system is still slightly 
out of equilibrium. However in the next section we shall see that a direct computation of 
the non equilibrium properties can be done. 



3 The non equilibrium equations 

It was found quite recently ||^, |^ that the non equilibrium behaviour of the system can be 
described directly for the infinite system (i.e. after having taken the limit oo), by 

introducing the average correlation function and the response function defined as 

Cit,t')= hm 



N^oo N 

G(M') = lim (14) 

In the equilibrium regime time translation invariance implies that these functions depend 
only on the time difference. We consider here the case where the system at time zero 
starts from a random configuration. Only positive times are possible and time translation 
is explicitely broken. 

Closed equation can be written for these two functions. They are 



dt' 
dG{t, t') 

dt' 



Eg[C,GI (15) 



where Ec[C.,G\ and Eg[C,G\ have an explicit form (non local in time) which depend on 
the problem. In some case one can expand Ec[C, G] and Eg[C, G] in powers of C and G. 
The solution of these equations can be computed numerically and one can obtain a great 
amount of information in this way. 

From the analytic point of view one can study these equation in the adiabatic approx- 
imation, where one set to zero the time derivative. This approximation is justified in the 
large time region. The resulting equations are non trivial; they are 

Ec[C,G]=Eg[C,G] = Q. (16) 
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The solution of these equations can be simphfied 0, H, P by noting that they are reparametriza- 
tion invariant, i.e. if C and G are a solution also the functions 



Ci,{t,t') = C{h{t),h{t')), 
G,it,t') = G{h{t),h{t'))^, (17) 

are an other solution of the adiabatic equations, for an arbitrary choice of the function h. 

Reparametrization invariance strongly simplifies the study of the adiabatic equations 
and many results can be obtained in this limit. In some case it can be proved that for large 
times the internal energy tends to the equilibrium value, while in other case finds that 
there is dynamical transition at a temperature Tq. At lower temperature one finds that 
the dynamical energy does not tend at large times at the equilibrium value and therefore 
metastable states are present. 

The evaluation of the reparametrization invariant quantities morally corresponds to the 
evaluation of the properties of the solutions of the mean fields equations of the previous 
section, although it contains more information. 

The more difficult part, which at the present moment we can do only numerically, 
consists in computing quantities that are not reparametrization invariant, as the time 
dependence of the energy. This computation morally corresponds to the evaluation of the 
barriers separating the solutions of the mean fields equations of the previous section and 
it is not a surprise that it turns out to be much more difficult. Technically one ends up 
with a well defined and difficult mathematical problem, very similar in spirit, but more 
complicated, of the non linear velocity selection problem, which has been widely studied 
in the past. 

It is extremely satisfactory that the very difficult problem of computing analytically 
the non equilibrium dynamics in these systems is now under control and I am convinced 
that the mathematical difficulties may be surmounted, may be with some help from our 
more mathematically minded friends. 

The results obtained from this dynamical approach have the advantage to be easily 
compared with those obtained experimentally in spin glasses, where the condition N very 
large with respect to t is certainly satisfied. A very interesting phenomenon which appears 
is aging, i.e. the dependence of the experimental results on the age of the system [^, |TT|, . 
A detailed discussion of this point would make this note too long. 



4 Glasses in the mean field approximation 

In spin glasses the Hamiltonian is random as an effect of quenched random disorder. In real 
glasses the Hamiltonian is not random and the quenched disorder is dynamically generated 
at low temperature. We can ask how much of the qualitative and quantitative results which 
have obtained in spin glasses may be transferred to glasses. 

In order to understand this point we have started to study models in which the Hamil- 
tonian does not contain quenched disorder and to compare the results with those coming 
of random Hamiltonian [|13| . 



Our strategy is the following. We want to study the properties of a given Hamiltonian 
Hg which is not random. We consider a class of Hamiltonians Hr, of which Hq is a 
particular case. We choose the class if/j in such a way that the statistical properties of 
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Hg and that of a generic Hamiltonian in Hji are as similar as possible. In the best case 
we can obtain that the two corresponding free energies coincide in the high temperature 
expansion 0. 

After having constructed Hr in an appropriate way, we can suppose that the statistical 
properties of Hq and are the same or, if they are different, we can construct a per- 
turbative expansion which compute this difference. It is clear that this approach may be 
successfully in the high temperature region (more or less by construction) and it may also 
reproduce the behaviour in the glassy region, included the dynamic and static transitions. 
However it is cannot certainly reproduce the possible existence of a crystal phase. 

I will present now a simple model in which this approach works very well at all the 
temperature and it misses only the crystal phase, which exists only for intentionally chosen 
Hamiltonians Hq |]14|- The Hamiltonian is is the same as in eq. (|I|). 



In the case of Hq we have 



J,. = iV-V2sin(^), (1^ 



while in the case of we have that J is a random orthogonal symmetric matrix, i.e. a 
random symmetric matrix which satisfies the constraint 

'^Ji,kJk,j = (19) 
k 

It is easy to check that the first Hamiltonian is a particular realization of the second one. 

A details computation (partially analytic and partially numeric) for the Hamiltonian 
Hg shows that in the low temperature phase there are two different limits when N goes to 
infinity, one for generic N and an other for N odd, such that p = 2N + 1 is prime. Only 
in this second case there is a crystalline phase at low temperature. At all temperatures 
the model for generic N behaves in the same way as the model described by and it 
undergoes a replica symmetry breaking transition. On the contrary the dynamical glassy 
transition is present for all N. The strategy of computing the properties of a given system 
by approximating it by a random system works very well in the glassy region and it misses 
the crystal phase, which exists only for non generic values of A^. 



There are many other system which can be studied using this stategy [|T5|, |I6|, but I 
will not discuss them. 



5 Conclusions 

We have seen that we begin to control the off equilibrium dynamics in many model systems. 
There are cases in which we have a glass transition with metastability. This feature is 
present only in the mean approximation, which is correct in the infinite range limit. 

A crucial and open problem is to understand how this feature of the infinite range 
approximation appears in more realistic finite range models, in which metastable states 
cannot exists. Let us describe a possible scenario for the (fragile) glass transition. 

The correlation time r diverges algebraically at Td, which is not the situation for real 
glasses where a divergence of the form 

r oc exp( ^ ) (20) 
(^ -JgJ 

^In general the behaviour of the system can be better controUed in the high temperature phase. 
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is observed. 

On the other hands metastabihty is present in the mean field approximation as soon as 
T < Tp. Generally speaking metastable states do decay in short range model by tunnelling 
effects. If we apply these ideas to the present case we find that 

r(xexp((7(T)), (21) 

where g(T) is a smooth function, which quite likely diverges at the temperature Tr as 

g{T)cx{T-Tnr''. (22) 

From this point of view the glass transition temperature should be identified with the 
temperature at which in the mean field approach the replica symmetry is broken. The 
dynamical transition found in the mean field approach is no more a real transition and 
denotes the onset of very slow dynamics. 

This line of thinking should be seriously investigated and one should develop the theo- 
retical techniques needed to compute the function g(T). The dependence of the exponent 
7 on the dimensions of the space is a rather interesting issue, which unfortunately it has 
not been studied up to now. Experiment on four dimensional glasses are notoriously rather 
difficult to be done, however the comparison of numerical simulations for three and four 
dimensional glass models should give rather interesting information and may be crucial 
to test further developments of the theory. For example the theory could predict that in 
(i-dimensions 7 = c? — 2, i.e. 

T oc exp(--- — ^ , , , ^, ). (23) 
(T - Tg)('^-2) ^ 

The only possibility for a new testing this kind of predictions (apart from the known three 
dimensional case) would be to do accurate numerical simulations in higher dimensions. 

The study of off equilibrium dynamics in glass related models has started only recently. 
Very interesting results have been obtained and I am confident that this approach will lead 
to a better understanding of real glasses. 
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